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Abstract 



The CDF preliminary analysis on polarized charmonium production at mod- 
erate transverse momentum, ~ 4 — 20 GeV, severely challenges the color 
octet model (COM), which predicts quarkonium to be transversely polarized 
with increasing pT- Based on this data, we analyze the compatibility of the 
Tevatron and the photoproduction at HERA in the context of the COM. Due 
to the uncertainty on the extraction of non-relativistic QCD (NRQCD) ma- 
trix elements and a lack of complete next-to-leading order calculations, one 
cannot completely rule out the COM. Nonetheless, both collider experiments 
seem to push the input matrix elements to opposite directions, and the puzzle 
of quarkonium polarization remains unsolved. 
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I. INTRODUCTION 



The simplest mechanism based on perturbative QCD to explain quarkonium production, 
the color singlet model (CSM) [|l|, is not able to describe charmonium hadroproduction. 
This model underestimates J/ip production, both in the central and forward |Q rapid- 
ity regions. The data show that the bound state of heavy-quark pair is produced not only 
in the color singlet configuration, but there are additional states contributing to the final 
colorless vector meson. Based on the non-relativistic QCD model (NRQCD) 0, quarko- 
nium production is understood as two-step phenomenon: cc pair production at perturbative 
level and the subsequent evolution to colorless vector meson through soft gluon emission at 
the non-perturbative domain. This argument is supported by the fact that the cc pair is 
produced at distance l/rriQ, itlq standing for heavy quark mass, which is much smaller than 
1/ Kqcd, the typical QCD scale for bound-state system. According to the color octet model 
(COM) formulation [^], a generic 5- wave quarkonium state is described by the Fock state 
decomposition, schematically given by 

I^q) = 0(l)|Qgp5f ^]) + 0{v)\QQ['pf]g) + 0{v')\QQ['Sf^]g) + 

0{v')\QQfsf]gg) + ... (1) 

where w is a typical velocity of the heavy-quark pair. We use the usual spectroscopic notation 
^■^"•"^Lj, and the color state is indicated by (1) for singlet and (8) for the octet. 

In the first approximation, the QQ system is produced in a color singlet state, which 
already has the quantum numbers of the physical quarkonium. The octet contributions are 
suppressed by powers of v and as- The latter is due to the extra soft gluon radiation needed 
to produce the correct color and/or quantum numbers. In principle, the state ^Pj^'' can 
produce xj mesons or evolve nonperturbatively to a vector meson. 

At the partonic level, the inclusive if) (generically denoting the charmonium J/ip and 
iIj{2S)) is given by 

da{a + 6^^ + X)=^ dan{a + cc[n] + X){Ot), (2) 

n 

where cc[n] stands for the quark-pair in the generic state n. We denote o"„ as the cross 
section for the short distance c-pair production, which can be calculated perturbatively. 
The matrix elements of the transition cc[n] —>■ ip, (O^), cannot be calculated in the usual 
perturbation theory. Fortunately, they are assumed to be universal, and can be extracted 
from experiments. 

Of course, in principle one could argue that the dominant long-distance matrix element 
should be {^Si^^'^), the cc state already with the correct quantum number and color of 
vector meson. However, as already stated, from the Fermilab experiments the CSM itself 
cannot explain the transverse momentum pt of the inclusive reaction pp ipX. The CSM 
differential cross section behaves like da/dpT ~ I/Pt; falling much faster than the data. 
The Pt dependence can be fixed combining the octet states '^Si^\ ^Sq, ^Pj, which is order of 
mlv^ according to NRQCD expansion. Particularly, the ^Si^^ is fundamental for explaining 
a harder p^ spectrum. If nature favors the vector S-wave octet state to evolve to a vector 
meson, there is a strong consequence on quarkonium production. Because the cc bound-state 
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is originated from gluon jet, the COM predicts quarkonium to be transverse polarized on 
the hmit Am1/p\ <^ 1 Indeed, it has been shown that in this hmit, the charmonium 



fragmentation function Dg^^ 
30(8) 



could be recovered 



The ^5r' plays the major role on the explanation of charmonium data at Tevatron, 
however the same is not true at HERA for z > 0.2, where z = ■ Pp/p.^ ■ Pp. In this 



kinematic region, the COM predictions for photopro duct ion [10-12] are dominated by the 



states ^5*0 and ^Pj. Moreover, in the Ref. |T^, the authors show that the polarization of 
J/ip produced from photoproduction depends only on the rate of the matrix elements of 
these states. 

Due to the universality, the values of matrix elements extracted from CDF data should 
reproduce HERA data. However, the COM predicts an excess of events compared with 
HERA data |p!3| , p!^ for z ^ 1. This discrepancy could be explained by the higher-order 
QCD effects [|15[ or by the intrinsic transverse momentum of the partons |T^. Nonetheless, 
it seems these tentative solutions cannot reproduce quite well not only z, but other relevant 
kinematic distributions 



In the following, we make a quantitative study of the COM in the light of the experimental 
data and we show that is quite difficult to accommodate the J/^p production at Tevatron 
and HERA, simultaneously. Even with the introduction of higher order QCD corrections, 
the COM will face another challenge: the interpretation of charmonium polarization. The 
preliminary CDF analysis [|l8l is pointing to unpolarized ip production, contradicting the 
COM predictions. 

Our strategy is at follows. We extract independently the non-perturbative matrix el- 
ements from both Tevatron and HERA data. After determining the solution that could 
satisfies both data, we show that actually it is incompatible to the polarization data. 



II. HADRO AND PHOTOPRODUCTION IN THE COLOR OCTET MODEL 



At the Tevatron, the inclusive production cross section can be written as the usual 
form. 



fa/p{Xa)fb/p{Xb) ^ab^ipxi^)^ 

where the a is given by Eq. (0) and at perturbative level 



(3) 



da. 



(A) 

a6— >ccX 



n 



dt 



Aab[n] + Bab[nM\) ■ A;,]^ + C,4n][e(A) ■ h]' + Dab[nM\) ■ A;J[e(A) ■ h], (4) 



where ka and kb are the momenta of the initial partons a and b and e(A) the polarization 



vector of ip. The complete analytic expression can be found in |p|,pJ]|. The sum over A yields 
the unpolarized cross section. 

Since we want to detect a vector meson, the lowest order at the partonic level should be 
2^2 process. At the Tevatron, the most important contribution comes from the subprocess 



g g cc[n\ g, 

although gq and qq bring some contribution, especially for high p^- 



(5) 
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At HERA, there are two types of mechanisms contributing to distinct regions of z. Direct 
i\} photoproduction, given by the partonic level subprocesses 



7* 5( cc\n\ g (6) 
7* cc[n] q{q), (7) 

which are important for z > 0.2, and by resolved photon mechanism through the partonic 
content of the 7*, dominant for small z. All the relevant analytic expressions are listed in 
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The complete ep ipX can be written as 

daep^ej/^xis) = r(Q^ y) da^^p^j/^xiW^), (8) 

where 

The dalW"^) can be related to the partonic cross section for the resolved photon process 
according to 

da^,p^j/^x{W^) = J fi/y{xi)fb/p{xb) d&ib^^xis)- (10) 
For i = 7, fif-yixi) = 5(1 — Xi) reproduces the expression for direct production. 

III. DISCUSSIONS AND CONCLUSIONS 

We have performed a numerical calculation of the charmonium production cross section 
fixing rric = 1.5 GeV and choosing the renormalization and factorization scale to be /i = 
Ptv "^V" '^^^'^6 "^V' ~ 2mc. We make our analysis for three different parton distribution 



functions (PDF's) in the proton; MRS (R2) 0, CTEQ 4L 0, and GRV 94 LO [p. For 



the resolved photons, we use GRV distribution function p2|. As the shape of cc[^Sq^'^] and 



ccfPj^^] are almost identical, we use the usual combination {^S^^^) + ^(^Pj^^) = Mk, fixing 
k = 3. 

Evidently a careful analysis on fi, as well as rric dependence would bring to better control 
of the theoretical uncertainties. Overall, even with an uncertainty with factor two, our 
conclusions still remain valid. Nevertheless, we consider a case with rric = 1-3 GeV for a 
more complete check. 

As we point out in the introduction, we extract the non-perturbative matrix elements 
independently, for both Tevatron and HERA. From the Fig. |l| we can see that the COM can 
accommodate quite well the CDF central {\ri^\ < 0.6) direct J/ip production data 0. We 
should emphasize that the same set of matrix elements brings to an extraordinary agreement 
with the DO forward (2.5 < \ri^\ < 3.7) J/ip production data |Q. However, this is not a 
surprise, once we only fit these free universal matrix elements without any constraint. 

In the Fig. |^ (|^) we show the z [y*, the rapidity of J/tp in the 7*p center-of-mass frame) 
distribution for J/tp production at HERA, and once again, it seems COM can in principle 
fit well the HI data IT 
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In the Table | we collect the best set of the color octet NRQCD matrix elements that fit 
both data set independently for the three PDF's we are considering. For the color singlet 
contribution, we have used = 1.2 GeV^, following 0. Since MRS (R2) is calculated 

at next to leading order, it is not surprising that it gives a different result compared to the 
leading order (LO) ones. Our numbers confirm the early results pointing out that at LO 
the COM has trouble explaining both data simultaneously. 

As we mention in the Introduction, this anomaly may be cured by adding corrections 
due to the intrinsic transverse momentum of the partons [|16| or higher-order (HO) QCD 
effects |T^. For small values of pt, the multiple-gluon radiations from the initial and the 
final state at the Tevatron become sizeable. In the [^] these corrections were estimated by 
Monte Carlo simulation using PYTHIA and the HO QCD could be parameterized as a 



K factor dependent on pt [|T5[- In fact, such corrections produce 

f^f)) = (0.47 ±0.09) X lO-^GeV^' 
M3 = (0.63 ± 0.34) X lO^^GeV^ 

for MRS (R2) PDF. This lower value of M3 brings to a better agreement with the one 
extracted from HERA experiments. 

In the Fig. ^ we display a parameter space for the color-octet NRQCD matrix elements. 
Although at 68% C.L. we still observe discrepancy between the bound for Tevatron HO 
QCD corrections and HERA, this picture changes dramatically at 95% C.L. We see that 
HERA favors much higher values for {^s[^^) than Tevatron, however they are not severely 
constrained. The reason is that the state ^S[^^ contributes only through resolved photon 
processes in a region where z < 0.4 , roughly speaking. In this intermediate region, the 
color-singlet contribution has a major role. Actually, this is clear if we consider a different 
c-quark mass value. For rric = 1.3 GeV, the state ^Si^^ has a bigger contribution, much 
closer to the experimental data. This means the CSM itself could describe the HERA data 
reasonable well, except in the region z ^ 1. 

The main conclusion from the Fig. |^ is that the introduction of HO QCD corrections 
brings a match between COM predictions at HERA and Tevatron, as already pointed in 
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With the extraction of the NRQCD matrix elements, we are now in the position to 
discuss the implication of these results on the charmonium polarization predicted by the 
COM. 

The quarkonium polarization can be measured from the angular dependence on ■?/' — *■ 

(xl+acos^e, (11) 



dcosO 

with a = (1 - 30/(1 + 0^ ^ = dcr^^V da\ 

From the expressions in [^, it is possible to calculate the quarkonium cross section for 
each polarization through the Eqs. (3) and (HI). Writing the polarization vector of quarko- 



nium in the recoil (s-channel helicity) frame EH], we found the pt dependence on a, as 



displaced in the Fig. ^. The polarizations were calculated for seven px bins, specified in [18 



At this point, we should be careful in comparing our results with the CDF preliminary 
analysis, once the data contain feed-down from Xc, which contributes to ~ 35%, and ip{2S) 
decay to J/ip. 
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In order to avoid these extra contributions to the prompt charmonium production, we 
also performed analysis on ip{2S) productions, which do not receive feed-down contributions, 
neither decay of the higher states. Although from the theoretical point of view the 'ip{2S) 
state is simpler to analyze, the available data are more limited statistically. 

We extracted the NRQCD matrix elements in a similar way we have done for J/tp. For 
the MRS (R2) PDF, we found (^sf ^) = (0.14 ± 0.03) x lO'^ GeV^ and very small value for 
Ms, compatible with zero. Following [|l^], we calculated the ip{2S) polarization for three pt 
bins, and once again the charmonium was found to be transverse polarized, as we can see 
from the Fig. ^ 

The HO QCD corrections, that worked well to solve the Tevatron/HERA discrepancy, 
actually worsen the already poor LO predictions, as the Fig. ^ indicates. 

In order to satisfy polarization data, the {^S[^^) contribution must vanish, since the 
M3 brings to a almost unpolarized ip production. Of course, there is a penalty doing just 
adjustment by hand. The nice fit, e.g., Fig. |l|is no longer held. Besides, the M3 value should 
be increased in order to have a better fit. From the Fig. ^, we see that this scenario will be 
disastrous if we compare with HERA bounds. 

Although there are strong evidences that COM is not working well, before ruling out 
this model, we should investigate the possible solutions to solve the paradox of charmonium 
production: 

• The complete QCD higher order corrections, which is not available so far, and con- 
tributions from higher cc states could in principle bring to a drastic change in the 
scenario. 

• The emitted gluons, in order to produce the physical quarkonium, are not so soft. 
Therefore, the polarization of the cc system is not conserved during the evolution to 
non-perturbative regime. 

On the other hand, although it is a strong statement, we could argue that the evolution 
of QQ system to a physical vector boson is not well understood; the splitting between per- 
turbative and non-perturbative regimes cannot be done trivially. This means that NRQCD 
is not appropriate to describe quarkonium production. Actually, if we remember that for the 
charmonium the perturbative expansion is based on 0{mcVc), maybe rric is not sufficiently 
small to allow such expansion. A closer examination on bottomonium states will be crucial 
to check if this state is held or not. 

Although is still early to make any strong conclusions, it seems that the COM is once 
again in trouble. At least in leading order cannot explain simultaneously the Tevatron and 
HERA data. The existing solution, the implementation of HO QCD corrections, worsen the 
strong prediction of this model, the production of transverse polarized quarkonium. 
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TABLES 



Jsf) 



MRS (R2) 



CTEQ 4L 



GRV 94 LO 



0.70 ± 0.17 (45 ± 29) 



0.54 ± 0.12 (25 ± 22) 



0.57 ± 0.12 (24 ± 22) 



+ 



4.85 ± 0.95 (0.39 ± 0.18) 



2.28 ± 0.55 (0.29 ± 0.14 ) 



2.07 ± 0.53 (0.30 ± 0.14) 



TABLE I. Leading-order Color-octet NRQCD matrix elements in units of 10 ^ GeV^ for direct 
J/ip production at Tevatron (HERA). 
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FIGURES 
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FIG. 1. The Pt distribution data (circles) for direct forward J/^p production < 0.6) 

from the CDF Collaboration at ^/s = 1.8 TeV. The solid curve represents the COM prediction 
after choosing the values for the NRQCD matrix elements given in Table | for the CTEQ 4L parton 
distribution function. The dashed curve shows the color singlet contribution. 
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o HERA HI 




FIG. 2. The z distribution for the inelastic J/V' production at HERA from HI Collaboration 
in the kinematic region 4 < < gg QeV^, V%^ji^ > 4 GeV^, AO < W < 180 GeV and z > 0.2. 
The solid curve represents the COM prediction after choosing appropriate values for the NRQCD 
matrix elements given in Table | for CTEQ 4L parton distribution function. The dashed curve 
shows the color singlet contribution. 
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o HERA HI 
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FIG. 3. Same as in Fig. || for the y* (the rapidity of J /if: in the 7*p center-of-mass frame) 
distribution. 
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FIG. 4. Parameter space for the Color Octet NRQCD matrix elements. The bounds on {^S^"^) 
and M3 for Tevatron and HERA are displaced at 68% C.L. (solid lines) and 95% C.L. (dashed 
lines). The results are for MRS (R2) parton distribution function. 
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FIG. 5. The polarization parameter a as a function of pT for the inclusive prompt J/ip pro- 
duction at the Tevatron. The bounds at LO (solid lines) and HO QCD (dashed lines) are based on 
68% C.L. including only errors due to the experimental data from CDF preliminary analysis for 
lyj/^l < 0.6. The results are for MRS (R2) parton distribution function. 
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o CDF preliminary 
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FIG. 6. Same as Fig. |5| for the ^^{23) production at the Tevatron. 
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